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To determine the Hilbert space and inner product for a quantum theory defined 
by a non-Hermitian PT-symmetric Hamiltonian H, it is necessary to construct a 
new time-independent observable operator called C. It has recently been shown 
' that for the cubic PT-symmetric Hamiltonian H = p 2 + x 2 + iex 3 one can obtain C 

perturbation expansion in powers of e. This paper considers the more difficult 
case of noncubic Hamiltonians of the form H = p 2 + x 2 (ix) s (5 > 0). For these 
Hamiltonians it is shown how to calculate C by using nonperturbative semiclassical 
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methods. 
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CT) ■ I. INTRODUCTION 

^ In 1998 it was discovered that the class of non-Hermitian Hamiltonians 

O 

g: H = p 2 + x 2 (ix) & (S>0) (1) 

has a positive real spectrum pj. It was conjectured in Ref. [IJ that the spectral positivity 
was associated with the space-time reflection symmetry (VT symmetry) of the Hamiltonian. 
The Hamiltonian (0) is VT symmetric because x — > —x and p — > — p under parity reflection 
V, and x — > x, p — > —p, and i —>■ — z under time reversal T. Many other "PT-symmetric 
quantum mechanical models have been examined and aproof of the positivity 

of the spectrum of H in (JTJ) was subsequently given by Dorey et al. (7(. 

Once the positivity and reality of the spectrum of a PT-symmetric Hamiltonian H, such 
as the non-Hermitian Hamiltonian (0), has been established one must then demonstrate 
that H defines a consistent unitary theory of quantum mechanics. To do so, one shows that 
the Hilbert space on which the Hamiltonian acts has an inner product associated with a 
positive norm and that the time evolution induced by such a Hamiltonian is unitary; that 
is, the norm is preserved in time |8j, |9j . Specifically, for a complex Hamiltonian having an 
unbroken VT symmetry, one must construct a linear operator C that commutes with both 
H and VT. One can then show that the inner product with respect to CVT conjugation 
satisfies the requirements for the theory defined by H to have a Hilbert space with a positive 
norm and to be a consistent unitary theory of quantum mechanics. [The term unbroken VT 
symmetry means that every eigenfunction of H is also an eigenfunction of the VT operator. 
This condition guarantees that the eigenvalues of H are real. The Hamiltonian in has 
an unbroken VT symmetry for all real 5 > 0.] 

In a conventional quantum theory the inner product is formulated with respect to ordi- 
nary Dirac Hermitian conjugation (complex conjugate and transpose). Unlike the situation 
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with conventional quantum theory, the inner product for a quantum theory defined by a 
non-Hermitian PT-symmetric Hamiltonian depends on the Hamiltonian itself and is thus 
determined dynamically. One must first find the eigenstates of H before knowing what the 
Hilbert space and the associated inner product of the theory are. The Hilbert space and 
inner product are then determined by these eigenstates. 

We emphasize that the key breakthrough in understanding these non-Hermitian VT- 
symmetric quantum theories was the discovery of the operator C p| . The problem is therefore 
to construct C for a given H. In Refs. @, Q it was shown how to express the C operator 
in coordinate space as a formal sum over appropriately normalized eigenfunctions <f) n (x) of 
the Hamiltonian H. These eigenfunctions satisfy the Schrodinger equation 

#0 n (x) = E n( j) n (x)i (2) 
and, without loss of generality, their overall phases are chosen so that 

VT(j) n (x) = <j> n (x). (3) 
With this choice of phase, the eigenfunctions are then normalized according to 

dx[<p n (x)f = (-iy\ (4) 

c 

The contour of integration C is described in detail in Ref. For the quantum theories 
described by the Hamiltonian ([TJ, C can be taken to lie along the real-x axis if 5 < 2. 

In terms of the eigenfunctions defined above, the statement of completeness for a theory 
described by a non-Hermitian PT-symmetric Hamiltonian is 

J2(-i) n Mx)My) = t(x-y) (5) 

n 

for real x and y. The unusual factor of (— l) n arises because of the convention (J3J). The 
coordinate-space representation of C is 



c 0> v) = 53 ( f ) n(x)(f)n(y)- (6) 

Only a non-Hermitian PT-symmetric Hamiltonian possesses a C operator distinct from 
the parity operator V. Indeed, if one evaluates the summation (JBJ for a PT-symmetric 
Hamiltonian that is also Hermitian, the result is V, which in coordinate space is d~(x + y). 

The coordinate-space formalism using (0) has been applied successfully to the cubic 
Hamiltonian 

H =\p 2 + \ii 2 x 2 + iex 3 , (7) 



and C was constructed perturbatively to order e 3 11 1. The approach in Ref. [ll| was to 
calculate the eigenfunctions as perturbation series in powers of e and then to substitute 
these series into (JBjl . The sum was then evaluated order by order in e. This techni que was 



also used to calculate C to order e for the cubic complex Henon-Heiles Hamiltonian |12j, [1 

H = \( P l + Pl) + \{x 2 + y 2 ) + iex 2 y, (8) 
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which has two degrees of freedom, and for the Hamiltonian 

H = \ (p 2 x + p 2 y + pl) + I (x 2 + y 2 + z 2 ) + iexyz, (9) 



which has three degrees of freedom [1J: 

Calculating the operator C by evaluating the sum in (jUJ) directly is difficult in quantum 
mechanics and hopeless in quantum field th eory because to do this it is necessary to deter- 
mine all the eigenfunctions of H. In Refs. 0, 3] we showed that there is a simpler way 
to calculate C based on three crucial properties of this operator. First, C is PT-symmetric 
(that is, it commutes with the space-time reflection operator VT), 

[C,VT] = 0, (10) 

although C does not commute with V or T separately. Second, the square of C is the identity, 

C 2 = 1, (11) 

which allows us to interpret C as a reflection operator. Third, C commutes with H, 

[C,H] = 0, (12) 

and thus is time independent. 

We also observed that there is a natural way to represent C as an exponential of a real 
Hermitian operator Q multiplying the parity operator V fl4l Ei|: 

C = e Q{x ' p) P, (13) 

where x and p are the dynamical operator variables. (This exponential representation was 
first noticed in Ref. ^]J.) The advantage of this representation is that when C is written in 
this form, properties (jlOJ) and (jllj) are equivalent to the symmetry conditions that Q(x,p) 
be an even function of x and an odd function of p. The remaining condition (fT2*|) yields an 
operator equation that determines Q(x,p). This operator equation can be solved perturba- 
tively for Hamiltonians having a cubic interaction term. The perturbative procedure is so 
easy that the C operator can even be found for cubic quantum field theories such as 3 3 

tf=I Vr 2 + I(V^) 2 + ±/iV + ^ 3 . 

To illustrate the perturbative procedure, we show how to find C for the Hamiltonian ((7|). 
We expand the operator Q(x,p) as a series in odd powers of e [3 Ell]: 

Q(x,p) = eQ^x.p) + e 3 Q 3 {x,p) + e 5 Q 5 {x,p) + ... . (14) 

We then substitute (fPfj) into (|T2*|) and collect the coefficients of like powers of e n for n = 
1,2,3, . . .. The result is a sequence of equations that can be solved systematically for the 
operator- valued functions Q n (x,p) {n = 1, 3, 5, . . .) subject to the symmetry constraints that 
ensure the conditions (|1U|) and (fTTj) . The first three of these equations read 

[H ,Qi] = —2Hi, 

[Ho,Qs] = -|[Qi,[Qi,#i]], 

[Ho,Qb] = ^[QuiQxAQuiQuHxM-KiQtAQs^Ht^ + lQs,^,^]]), (15) 
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where H = \p 2 + \p 2 x 2 and Hi = ix 3 . We now substitute the most general polynomial 
form for Q n using arbitrary coefficients and then solve for these coefficients. For example, 
to solve the first of the equations in (|15|). we take as an ansatz for Qi the most general real 
Hermitian cubic polynomial that is even in x and odd in p: 

Qi(x,p) = Mp 3 + Nxpx, (16) 

where M and N are undetermined coefficients. The operator equation for Qi is satisfied if 

M = -f/x- 4 and N = -2p~ 2 . (17) 

In Ref. we determined the C operator to seventh order in perturbation theory. We were 
able to perform this high-order calculation because our procedure does not make use of the 
summation in (pj). and therefore we did not need to find the wave functions <p n (x). 

Unfortunately, this perturbative procedure for finding the operator C is ineffective when 
the Hamiltonian is not cubic because this approach leads to complicated and unwieldy 
infinite sums |l4| . Furthermore, the perturbative procedure does not work in the massless 
limit. [Note that the coefficients in (JT7j) diverge when p = 0.] Thus, a more powerful method 
is needed to find the C operator associated with the Hamiltonian 

In this paper we introduce a completely new approach based on nonperturbative semi- 
classical methods rather than perturbative methods. The procedure is straightforward; the 
WKB physical-optics approximation to the eigenfunctions <f> n (x) is substituted into the for- 
mula © for the operator C and the summation is performed. This procedure is justified 
because the WKB approximation is asymptotically accurate in the limit of large n. We 
will see that when the WKB eigenfunctions are used to evaluate this sum, the result is a 
singular operator. The error incurred by including the small-n terms in the sum is finite, 
and it therefore may be neglected. 

In the next section we show how to perform the WKB calculation of the C operator and 
in Sec. IHII we discuss the properties of our solution for C and consider the possible extension 
of this work to quantum field theory. 

II. WKB CALCULATION OF THE C OPERATOR 

We begin this section with a review of the WKB analysis of the conventional two-turning- 
point problem. We then generalize the standard treatment of the two-turning-point problem 
from the real-x axis to the complex plane, where we apply it to the Hamiltonian JT]). 

A. WKB on the Real Axis 

The conventional treatment of the two-turning-point eigenvalue problem on the real axis 
is described in Ref. jl6j . This problem is expressed in terms of the differential equation 

- e 2 0" (x) + V(x)<f>(x) = E<f>{x), (18) 

where e < 1 is a small positive perturbation parameter. The parameter e is used to or- 
ganize the WKB expansion and to distinguish between different orders in the semiclassical 
perturbation theory. We assume that the potential V(x) rises as x — >• ±oo and that there 
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are two turning points at x = A and at x = B with A < B. The turning points A and B are 
solutions to the algebraic equation V(x) = E. The classically allowed region is A < x < B 
and the classically forbidden regions are x < A and x > B. The quantization condition on 
the eigenvalue E in (|18|) is the requirement that the wave function <p{x) vanish as x — > ±00. 
This quantization condition leads to a discrete spectrum E n (n = 0, 1, 2, • • • ). 

For fixed n and small e, WKB theory gives a good approximation to the eigenvalues E n 
and the corresponding eigenfunctions (f) n (x). (For fixed e the WKB approximation becomes 
accurate as n — > 00.) To leading order in WKB theory (the physical- optics approximation) 
the quantization condition for the eigenvalues E n reads 

dx ^E n - V(x) ~ (n + ~J vr (e -> + ). (19) 

Also, the physical-optics approximation to the wave function 4> n (x) in the classically allowed 
region A n < x < B n is given by 



(f> n (x)~C[E n -V(x)}- L/ *sm 

The complete WKB physical-optics approximation to the wave function <f) n (x) actually con- 
sists of five separate formulas, but to calculate the operator C we will not need the approx- 
imations to the wave function in the regions near the turning points at A n and B n and in 
the two classically forbidden regions. (These formulas are given in Ref. [16].) 

The multiplicative constant C in (|2(J|) is determined by the usual normalization condition 
that the integral of the square of the eigenfunction be unity: 

/oo 1 pB 

dx<j)l{x) ~ -C 2 / dx [E n - V(x)]- 1/2 (e - 0+). (21) 
-00 * J A 

The calculation that leads to this asymptotic evaluation of the integral is explained in detail 
in Ref. 0] . Note that the condition in (}2~T]) employs the conventional normalization criterion 
that is used for Hermitian Hamiltonians rather than the normalization used in (jlj). 




- ^ dty/E n -V(t) + j 
e JA„. 4 



( e ^o + ). 



(20) 



B. WKB in the Complex Plane 

The Schrodinger eigenvalue equation (0) associated with the Hamiltonian is 

- e 2 £(x) + x 2 (za;)Vn(x) = E n <f> n (x). (22) 

We have inserted the small positive parameter e to organize the structure of the semiclassical 
approximation but at the end of the calculation we will set e = 1 and use the fact that WKB 
becomes accurate as n — » 00 to justify our results. 

To construct a physical-optics approximation to the eigenfunctions 4> n (x) for this equa- 
tion, we begin by finding the turning points A n and B n , which satisfy E n = x 2 (ix) s : 

A n = E^e' in+ ^ and B n = E^ e~$* . (23) 

These turning points lie on the real axis when 5 = and rotate downward into the complex-x 
plane as 5 increases from 0. (The turning points are shown in Fig. [I]) 
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The WKB analysis is done along a contour passing through these turning points. This 
contour, which is shown as a solid line in Fig. ^ lies in the lower-half complex-x plane in 
the asymptotic Stokes wedges in which the eigenf unctions vanish exponentially as \x\ — > oo. 
For large |x| the WKB contour is asymptotic to the centers of the wedges; these centers are 
shown on Fig. ^as dashed lines. In the quadrant Hex > 0, Imx < the center of the wedge 
lies at the angle — g^ 71 "; the upper and lower edges of the wedge (not shown in Fig. |TJ) 
lie at the angles — ■^p^n and — ^d^tt. A detailed description of the WKB contour and the 
asymptotic wedges is given in Ref. 

The WKB formulas ()19|) - (|21|) for the eigenfunctions 4> n {x) are valid in the comp lex- 
re plane as well as on the real-x axis. Our ultimate objective is to use these formulas to 
evaluate the sum (jHI) and thereby to determine the C operator. However, to demonstrate the 
techniques needed to evaluate this sum we will first show how to evaluate the sum in (jSJ) that 
expresses the completeness condition and we will verify the (approximate) completeness of 
the WKB wave functions in (J2U|) . 

We begin our analysis by finding the WKB approximation to the product <j) n (x)<j) n (y). 
We will then perform the sum in (jSJ) under the assumption that the arguments x and y of the 
eigenfunctions are real and of order 1. We must take x and y to be real because, as we will 















x < 


] y 







FIG. 1: Schematic representation of the contour (solid line) in the complex- x plane along which 
the WKB approximation to the Schrodinger equation (|22|) is derived. This contour passes through 
the two turning points A and B given in (J23|). These turning points are indicated by dots. For 
large |x| the WKB contour approaches the centers of the asymptotic wedges (dashed lines) in 
which the eigenfunctions vanish exponentially as \x\ — > oo. The eigenfunctions are oscillatory on 
the WKB contour between A and B. (This region is the complex generalization of the classically 
allowed region in conventional WKB theory.) Outside of this oscillatory region the eigenfunctions 
decay exponentially along the WKB contour. (These regions are the complex generalizations of 
the classically forbidden regions.) The eigenfunctions are analytic in the complex- x plane except 
on the branch cut, which runs along the imaginary- x axis from the origin to ioo. To calculate the 
C operator the WKB approximation (|20|) to the eigenfunction <j) n is evaluated on the real-x axis at 
the points x and y, which are at a distance of order 1 from the origin. 
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see, they will appear as arguments of a Dirac delta function, and the delta function is only 
defined for real argument. (The points x and y are shown in Fig. 0) We will assume that 
n>l, and therefore that E n 3> 1. Under this assumption we can make the approximations 



[E ri 



x 2 (ix) s 



-1/4 



£-1/4 



0" 



and 



where 



I dt ^/E n -t 2 (tty 



In + 



dt ^E n - t 2 (it) 1 



7T 



V ^+ 0/ g 4+2,5 _ 



Tp 4+26 

9 r f8+3<5\ 



(24) 



(25) 



(26) 



Using the trigonometric identity 2 sin (a + |) sin (& + f ) = cos(a — 6) + sin(a + 6), we obtain 
the WKB approximation to the product (fi n (x)(fi n (y): 



(j) n {x)(j)n{y) 



-l) n 



E n j2 n dt [E n -t 2 (ix) 



-1/2 



+ sin 



cos 



-\/E n {x + y) 



-y/E n (x 



cr 



y) 



(27) 



where the factor of (— l) n comes from the normalization condition (0J). 

We will see that this factor of (— l) n plays a crucial role in our evaluation of sums because 
oscillatory terms are suppressed relative to nonoscillatory terms. This factor arises because 
the phase-fixing condition © requires that the eigenfunctions <p n (x) be VT symmetric. The 
harmonic oscillator (5 = 0) illustrates clearly the origin of the (— l) n factor. Requiring that 
the harmonic oscillator eigenfunctions be VT symmetric (that is, that they be functions of 
ix) means that they have the general form (j) n {x) = i n e~ x ^ 2 He n (x), where He n (x) is the 
usual Hermite polynomial. It is the extra factor of i n that gives rise to the factor of (— l) n 
in the normalization condition The factor (— l) ra persists for all 5 > 

Next, we use the quantization condition (fTT?|) to obtain the identity 



n + 



1 



1 



7T 



(In 



(28) 



where * denotes complex conjugation. This identity allows us to rewrite (j2*7|) as 

1 



<j>n(x)<j> n {y) 



E n f* n dt \E n -t\ix) 



-1/2 



(-l) n cos 



+ COS 



E n (x + y) 



E n (x - y) 



0" 



(29) 



The next step is to perform the sum over n. We do so by converting the sum to an 
integral using the general WKB density of states result 



1 f Bn 

dE n — dx [E. 



V{x)\ 



-1/2 



tt dn 



(30) 
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which is obtained by differentiating the quantization condition (|19|). Note that one must 
differentiate the endpoints as well as the integrand with respect to E n , but the endpoints 
give a vanishing contribution because E n — V(A n ) = E n — V(B n ) = 0, where the potential 
V{x) = x 2 (ix) s . 

Let us substitute (|HU|) into (0) and examine the sum over n of the first term in the curly 
brackets in (|2T?|) . Letting r = y/E^/e, the sum over n in (JHJ) becomes the integral 

1 f°° 

— \ dr cos [r{x — y)] = S(x — y), (31) 

7T Jo 

which is precisely the statement of completeness. Of course, this completeness result is only 
an approximation because the WKB approximation is only accurate for large n. Thus, the 
early terms in the sum over n introduce an error in this result. This error is equivalent 
to changing the lower limit in the integral in (|31|) from to some other finite number L. 
However, the error that arises when is replaced by L is finite, and when x is near y, this 
error is small and may be neglected in comparison with 5(x — y), which is singular at x — y. 

The second term in the curly brackets in ()29|) does not contribute to the sum over n in (JSJ) 
because there is a factor of (— l) ra in this term. The rapid oscillation produced by this factor 
causes a cancellation between successive terms in the summation. As a result, this term 
contributes negligibly to the sum and the delta function result in (|31|) remains unchanged. 

Next, we turn to the evaluation of the C operator in (jUJ). Because there is no factor of 
(— l) n in this sum, it is now the first term in curly brackets in ()29|) that is unimportant and 
it is the second term that contributes. Converting this sum to an integral using (}3T)j) . we 
obtain in coordinate space 

C(x,y) = - / dr cos r(x + y) + iKr^e^ , (32) 
7T Jo L J 



where we have used the result that 



4+<5 

I n -r n = iKEn +2S (33) 



with the n-independent constant K given by 



-^i4^1. (34) 



For small e the cosine function in ()32j) can be expanded into a two-term Taylor series and 
the first integral can be done as in (J31j) : 

C(x, y) = 5(x + y) — iKe~ 6 — \ dr r~ & sin [r(x + y)] . (35) 

^ Jo 

To evaluate the remaining integral in (J33j) . we transform to momentum space: 

C(p,q) = j ' dxe~ ipx j ' dye iqy C(x,y). (36) 
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The Fourier transform of the first term in ()35|) is 2n5(p + q), which is just the parity operator 
V(p, q) in momentum space. The Fourier transform of the second term in (|35|) is 



—iKe 2 + s dxe 



-ipx 



dye 



iqy 



7T 



2z 



-^e^ y dx j dy J drr^l ^i(r-P>+i(<l+r)y _ e 
9(p)p^ - 6{— p)(-p)^ 5(p + q), 



-ir(x+y) j 
-i(p+r)x+i(q— 



(37) 



where 0(t) = (t < 0), 0(t) = 1 (t > 0) is the step function. Thus, our result for the WKB 
approximation to the C operator in momentum space is 



C{P, q) = r[p, q) — e 2 + 5 _ /oin J s sin 



1 



4 + 25 



4+5 4+0 

9(p)p— - 6(-p)(-p)*+* V(p,q). (38) 



4+6 " 



Finally, we rewrite this result in the exponential form C = e®V in (j!3j) . set e = 1, and 
identify the function Q: 



1 U+2<W 



7T 



5 



sin 



4 + 25 



4+5 



(39) 



This is the principal result of our investigation. 



III. DISCUSSION AND CONCLUSIONS 

We have found a physical-optics WKB approximation to the C operator in momentum 
space for the Hamiltonian in Note that the results in ()38)1 and (}39|) are valid for all 
5 > 0. This represents a major advance over previous work; previously, the C operator 
was only determined for the special case of a cubic interaction term 5 = 1. The methods 
used earlier were perturbative and were not as powerful as the nonperturbative semiclassical 
methods used in this paper. 

The function Q(p, q) exhibits the required symmetry properties mentioned earlier, namely, 
that it be an even function of x and an odd function of p. Furthermore, because the sine 
function in (|3*5|l and (|3*9*j) vanishes when 5 = 0, we see that the C operator becomes identical 
with the parity operator V when the Hamiltonian is Hermitian in the conventional sense. 
The vanishing of Q for the 5 = case can be traced back to the fact that for the harmonic 
oscillator the integral I in (JUJ) is real. Thus, I - I* = in ((23) and 

Note that the results in (jSHJ) and (|3l?j) simplify substantially when 5 is an odd integer 
because the step function 9 drops out of these formulas. For 5 = 2k + 1 we have 

^ v^r(|±f) /7r(l + 2A;)\ 5 +2* 

Gfo *) = -T(ifr sm {^TW ) p3+2K (40) 

When k = (the case of a cubic theory) we have 

Q( P a) = _ ^(1/3) 5/3 
wq) 5r(5/6) P ' 
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and when k = oo we have the especially simple result that 

Q(p,q) = -2p. 

Another special case is 5 = 2. For this case the Hamiltonian ((TJ is quartic and 

We hope to generalize the advance reported in this paper to noncubic "PT-symmetric 
quantum field theories, such as a —g(p 4 theory. A — g<p 4 quantum field theory in four- 
dimensional space-time is a remarkable model because it has a positive spectrum, is renor- 
malizable, is asymptotically free [T3|, and has a nonzero one-point Green's function G\ = ((f). 
If we can perform the corresponding semiclassical approximation for a — g(p 4 quantum field 
theory, we may be able to elucidate the dynamics of the scalar sector in a variant of the 
standard electroweak model and possibly even predict the mass of the Higgs particle. 
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